Alternative interpretation of the sign reversal of secondary Bjerknes force acting 

between two pulsating gas bubbles 
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It is known that in a certain case, the secondary Bjerknes force, which is a radiation force acting 
between pulsating bubbles, changes, e.g., from attraction to repulsion as the bubbles approach 
each other. In this paper, a theoretical discussion of this phenomenon for two spherical bubbles is 
described. The present theory based on analysis of the transition frequencies of interacting bubbles 
[M. Ida, Phys. Lett. A 297, 210 (2002)] provides an interpretation, different from previous ones 
(e.g., by Doinikov and Zavtrak [Phys. Fluids 7, 1923 (1995)]), of the phenomenon. It is shown, 
for example, that the reversal that occurs when one bubble is smaller and another is larger than a 
resonance size is due to the second-highest transition frequency of the smaller bubble, which cannot 
be obtained using traditional natural-frequency analysis. 
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I. INTRODUCTION 

It is known that two gas bubbles pulsating in an acous- 
tic field undergo an interaction force called the secondary 
Bjerknes force This force is attractive when the 

bubbles pulsate in-phase with each other, while it is re- 
pulsive otherwise; that is, the phase property of the bub- 
bles plays an important role in determining the sign of the 
force. In a seminal paper published in 1984 |4|, Zabolot- 
skaya, using a linear coupled oscillator model, showed 
theoretically that in a certain case, the sign of the force 
may change as the bubbles come closer to one another. 
This theoretical prediction was ensured by recent exper- 
iments that captured the stable, periodic translational 
motion of two coupled bubbles [5J, resulting from the 
sign reversal of the force at a certain distance between the 
bubbles. Zabolotskaya assumed that this sign reversal is 
due to variation in the natural frequencies of the inter- 
acting bubbles, which results in shifts of their pulsation 
phases. The theoretical formula Zabolotskaya derived to 
evaluate the natural frequencies of two interacting bub- 
bles, which corresponds to one given previously by Shima 
@, is represented as 
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where Rio and R20 are the equilibrium radii of the bub- 
bles, loio and LO20 are their partial natural (angular) fre- 
quencies, u) is the angular frequency of an external sound, 
and D is the distance between the centers of the bub- 
bles. This equation predicts the existence of two natural 
frequencies per bubble, and is symmetric; namely, it ex- 
changes 10 and 20 in the subscripts of the variables to 
reproduce the same equation, meaning that the two bub- 
bles have the same natural frequencies. 

During the last decade, a number of studies regard- 
ing the sign reversal of the force have been performed 
JaOJE II 13 El 13 13 III Among them, Refs. |1 
llOl [jjj also considered the relevance of the change in 



the natural frequencies (or resonance frequencies |15| ) to 
the sign reversal. In the present paper, we focus our 
attention on this case, although it has been shown that 
other factors, such as the nonlinearity in bubble pulsation 
nil H2l Ha. Il4| and the higher-order terms appearing 
in the time-averaged interaction force |5J which has been 
neglected in previous works, can also cause the sign re- 
versal. 

In 1995, Doinikov and Zavtrak |9j, using a linear math- 
ematical model in which the multiple scattering of sound 
between bubbles is taken into account more rigorously, 
predicted again the sign reversal. They also asserted that 
this reversal is due to the change in the natural frequen- 
cies. They assumed that the natural frequencies of both 
bubbles increase as the bubbles approach each other, re- 
sulting sometimes in the sign reversal. When, for exam- 
ple, both bubbles are larger than the resonance size (i.e., 
the case of < u> and W20 < and the distance be- 
tween them is large enough, they pulsate in-phase with 
each other. As the bubbles approach each other, the nat- 
ural frequency of a smaller bubble may first, at a certain 
distance, rise above the driving frequency, and in turn 
the bubbles' pulsations become antiphase; the force then 
changes from attractive to repulsive. When, on the other 
hand, one bubble is larger and the other is smaller than 
the resonance size (e.g., lu±o > u> > CJ20) and the dis- 
tance between them is large, they pulsate out-of-phase 
with each other and the force is repulsive. As the dis- 
tance between the bubbles becomes smaller, the natural 
frequencies of both bubbles may rise, and when the nat- 
ural frequency of a larger bubble rises above the driving 
frequency, the repulsive force may turn into attraction. 
This interpretation is supported even in more recent pa- 
pers (3Gj 

Although this interpretation seems to explain the sign 
reversal well, it is opposed to the prediction given by 
Eq. m which reveals that the higher natural frequency 
(converging to the partial natural frequency of a smaller 
bubble for D — > 00 @, 0]) increases but the lower one 
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(converging to the partial natural frequency of a larger 
bubble for D — ► oo) decreases as the bubbles approach 
each other. 

In 2001, Harkin et al. |13| performed an extensive 
theoretical study concerning the translational motion of 
two acoustically coupled gas bubbles in a weak and a 
moderate driving sound field. Their theoretical model 
derived from first principles supports the experimental 
results by Barbat et al. §. In Sec. 7 of that paper, 
Harkin et al. also considered the influence of the change 
in natural frequencies on the sign of the force in order to 
explain the sign reversal for wio < u) and a->20 < Their 
explanation based on a formula given directly by Eq. Q 
is essentially the same as those by Zabolotskaya [4| and 
by Doinikov & Zavtrak HG3- 

The authors should note here that all the previous 
theoretical models mentioned above can describe (or ex- 
plain) the sign reversal. However, the interpretation we 
will provide in the present paper is different from the 
previous ones. 

The aim of this paper is to give an alternative inter- 
pretation of the sign reversal, one that may be more 
accurate than the previous ones that are based on the 
natural- frequency analysis. Recently, having reexamined 
the linear coupled oscillator model used frequently to an- 
alyze the dynamics of acoustically coupled bubbles (see 
Ref. [16j and references therein), we found that a bubble 
interacting with a neighboring bubble has three "tran- 
sition frequencies" , defined as the driving frequencies for 
which the phase difference between an external sound and 
the bubble's pulsation becomes ir/2 (or S nJ 2) , two of 



n/2), 



which correspond to the natural frequencies [16( . Among 
the three transition frequencies, the lowest one decreases 
and the remaining two increase as the bubbles approach 
each other. Meanwhile, for D — > oo only one of them con- 
verges to the partial natural frequency of the correspond- 
ing bubble. Namely, the transition frequencies defined 
as above are asymmetric. The use of the transition fre- 
quencies would allow us an accurate understanding of the 
sign reversal, because observing these frequencies pro- 
vides more detailed insights of the bubbles' phase prop- 
erties rather than that provided by the natural-frequency 
analysis. Using the theory for the transition frequencies, 
we arrive at a novel interpretation of the sign reversal. 



II. THEORIES 

In this section, we briefly review the previously ex- 
pounded theories regarding the natural frequencies, the 
transition frequencies, and the secondary Bjerknes force. 



can be represented as Rj = Rjo + ej(t) and \ej \ <C Rjo, 
where Rjo and ej are the equilibrium radius and the de- 
viation of the radius, respectively, and j = 1,2, ... ,N . 
The radius deviation can be determined by solving the 
linear oscillator model (see, e.g., Ref. [l7|). 
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where 
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are the partial natural (angular) frequencies of bubble 
j, Sj is the damping coefficient determined based on the 
damping characteristics of the bubbles 01 > Pdj is the 
driving pressure acting on bubble j, Kj is the effective 
polytropic exponent of the gas inside the bubbles, Pq is 
the static pressure, a is the surface tension, p is the den- 
sity of the liquid surrounding the bubbles, and the over- 
dots denote the time derivation. The driving pressure 
is represented by the sum of p cx and the sound pressure 
scattered by the surrounding bubbles, p s , as 

N 

Pdj = Pcx + Pa,jk- 
k=l,k=ij 

The value of p s j & is determined by integrating the 
momentum equation for linear sound waves, dp/dr — 
~pdu/dt, coupled with the divergence-free condition, 
d(r 2 u)/dr = 0, where r is the radial coordinate mea- 
sured from the center of a bubble and u is the velocity 
along r. Resultantly, the driving pressure is determined 
as 
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where Dj & is the distance between the centers of bubbles 
j and k. 

In a single-bubble case (i.e., for N = 1), Eq. © is 
reduced to 
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Assuming that p cx is written in the form of p ox = 
—P a sm.ujt (P Q is a positive constant), the harmonic 
steady-state solution of Eq. (QJ is given by 

ei = K S1 sin(wi - <j> S i), 



with 



A. Natural frequencies and transition frequencies 

Let us consider the linear volume oscillation of N— 
bubble system immersed in an incompressible liquid. 
Suppose that the time-dependent radius of bubble j, Rj, 



K S i = 
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From this result, one knows that the phase difference of 
<f>si = tt/2 appears (or, roughly speaking, the phase re- 
versal takes place) only at the natural frequency uj\q 0] , 
and the resonance response occurs at (or, more correctly, 
near) the same driving frequency. 
For N = 2, Eq. © is reduced to 
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where D = Di 2 = D 2 i- It is known that for a weak forc- 
ing (i.e., P a <§; Pq), this system has third-order accuracy 
with respect to 1/D, although it has terms of up to first 
order (the last terms) The harmonic steady-state 

solution for ei is 



where 



e\ = Ki sin(o>i — 4>i), 



0! - tan" 1 (^j G [0,2tt], 



As was proven in Ref. ^(|, this equation predicts the 
existence of up to three transition frequencies per bub- 
ble. Furthermore, as pointed out in the same article, the 
terms in the second (• • •) of Eq. © are the same as those 
on the left-hand side of Eq. JQ). These results mean that 
in a double-bubble case the phase reversal of a bubble's 
pulsation can take place not only at its natural frequen- 
cies but also at one other frequency. Because Hi ^ H 2 , 
Eq. (|SJ (and also Eq. J7J) is asymmetric, meaning that 
the bubbles have different transition frequencies. 

A preliminary discussion for a TV-bubble system ppj 
showed that a bubble in the system has up to 2N — 1 
transition frequencies, N ones of which correspond to the 
natural frequency. Namely, a bubble has an odd number 
of transition frequencies. This result can be understood 
as follows: Even in a multibubble bubble's pul- 

sation may be in-phase or out-of-phase with a driving 
sound 2 1 J when the driving frequency is much lower or 
much higher, respectively, than its natural frequencies; 
thus, in order to interpolate these two extremes consis- 
tently an odd number of phase reversals are necessary 

M. 



B. Secondary Bjerknes force 
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Mi = 8\w, M 2 = S 2 uj. 



Exchanging 1 and 2 (or 10 and 20) in the subscripts of 
these equations yields the expressions for bubble 2. 

The formula for the natural frequency, Eq. QJ, is de- 
rived so that Ki — ► 00 for 5i — > and 5 2 — > 0. Namely, 

f = LiL2 _^3b^ = . 

As mentioned already, this equation predicts the exis- 
tence of up to two natural frequencies in a double-bubble 
system. 

The transition frequencies of bubble 1 are determined 
so that 4>i becomes ir/2 (or 3n/2). Because F 2 + G 2 ^ 
[l6|. the resulting formula for deriving the transition 
frequencies of bubble 1 is 



H±F + M 2 G = Q. 
Assuming Si — > and S 2 — ► reduces this to 



(7) 



HiF = U 2 + ^ 2 



RioRw 41 ,, (V1 
LiL 2 j^— uj ) = 0. (8) 



The secondary Bjerknes force acting between the 
bubbles for sufficiently weak forcing is expressed with 

n nan, urn 
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l r 2 - ril' 



( 9 ) 

where Vj and rj are the volume and the position, re- 
spectively, of bubble j, (■ ■ •} denotes the time average, 
and || r 2 — r x \\ = D. The sign reversal of this force oc- 
curs only when the sign of cos(0i — 4> 2 ) (or of (ViV 2 ^j) 

changes, because Ki > and K 2 > 0. If the phase shifts 
resulting from the radiative interaction between bubbles 
are neglected, this force is repulsive when to stays be- 
tween wio and lu 2 q , and is attractive otherwise [lj . In the 
case where the radiative interaction is taken into consid- 
eration, the frequency within which the force is repulsive 
shifts toward a higher range, see, e.g., Refs. [9HlC|. 



The formulae reviewed above, except for that regarding 
the transition frequencies (Eqs. J7J and JSJ), are classical, 
and almost the same ones have previously been used in 
Ref. As will be shown in the next section, however, 
the following investigation based on Eq. J7J) coupled with 
Eq. @ gives a different interpretation of the sign reversal 
from the previous ones described using only the natural 
frequencies. 
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FIG. 1: Transition frequencies ui\ (rad/s) and W2 (rad/s) for 
Rio ~ 2 mm, ifeo = 5 mm, and the reduced damping, nor- 
malized by cjio (rad/s). The dashed lines show the transition 
frequencies that do not cause the resonance. 



III. RESULTS AND DISCUSSION 

In this section, we investigate the relationship between 
the transition frequencies and the sign of the secondary 
Bjcrknes force by using some examples. The first exam- 
ple is the case of Riq = 2 mm an d Ro n = 5 mm, which 
corresponds to a case used in Ref. [T(|. We assume that 
the bubbles are filled with a gas having a specific heat 
ratio of 7 = 1.4, and the surrounding material is water 
(a = 0.0728 N/m, p = 1000 kg/m 3 , P = 1 atm, and 
the speed of sound c = 1500 m/s). For the damping co- 
efficient, we adopt that used for radiation and thermal 
losses: 



w 2 Rjo 



•A 



(10) 



where the thermal damping coefficient /3th j and the effec- 
tive polytropic exponent determined by 0, 0, ^ 



Pth.j — a th,jj 
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Kj =7 

with 



(l + 4,i) H 



3(7- l)(sinhA-sinA) 
A(coshA-cosX) 



dth,j =3(7-1) 

A(sinh A + sinX) - 2(coshA - cosX) 



A 2 (coshA-cosA) 
X = R j0 (2u/ X G) 1/2 , 



•3(7- l)A(sinhA-sinX)' 



where we set \G = 2 x 10~ 5 m 2 s _1 . 

In order to clarify the following discussion, we first 
present results for the idealized condition of Sj « by 
resetting Sj — > Sj/100, and subsequently provide results 
given by the direct use of Eq. (|lf)|> . Figure ^ shows the 
transition frequencies of the bubbles, u>\ and u>2, calcu- 
lated using Eq. J7J) with the reduced damping, normalized 
by ujiq (= wi for 13 — ► 00). In those figures, I denotes 
the normalized distance defined as I = D/(R 10 + i?2o)- 
As mentioned previously, we can observe three transition 
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FIG. 2: fa/ir, fa/ir, and cos(<^i — (/>a) for (a) = l.Olwio 
(rad/s) and (b) u) = 1.03oj20 (rad/s), as functions of I. The 
dashed curves and the solid lines show the results given using 
the reduced damping and the real damping, respectively. 



frequencies, only one of which converges to lujq of the 
corresponding bubble for I — > 00. The second- highest 
transition frequency of bubble 2 is almost equal to the 
highest one of bubble 1; thus, the highest one of bubble 
2 is higher than that of bubble 1. The second highest 
one of bubble 1 and the highest one of bubble 2 do not 
cause the resonance response |16|. 

The dashed curves displayed in Fig. |2fa) show fa , (j> 2 , 
and cos(0i — fa), respectively, as functions of I. Here 
the driving frequency is assumed to be lo = l.Olwio, i.e., 
slightly above loiq. (In the present study, the driving fre- 
quency is set as w ~ u>io or lo ps LO20 so that the sign re- 
versal takes place at a sufficiently large I where the accu- 
racy of Eqs. JSJ and © is guaranteed.) As mentioned in 
Sec.[U it is known already that the sign reversal can take 
place when lo > u>iq > LO20 or loiq > lo > LO20 ; the present 
setting corresponds to the former case. We can observe 
one and two sharp shifts of fa and fa, respectively. At 
I « 3, both fa and fa shift almost simultaneously, but 
the sign reversal does not occur because the phase dif- 
ference fa —fa is hardly changed. At I « 13, only fa 
shifts, resulting in the sign reversal. In the former case, 
the phase shifts are caused by the natural frequencies. 
As mentioned previously, when Sj w both the bubbles 
have (almost) the same natural frequencies. The simul- 
taneous phase shift, thus, appears. The change of fa in 
the later case is apparently due to the highest transition 
frequency of bubble 2, which cannot be obtained by the 
traditional natural-frequency analysis. Namely, this sign 
reversal cannot be interpreted by using only the natural 
frequencies. 

We should note here that, to compute the phase delays 
fa and fa, we used the "atan2(o,6)" function in the C 
language, which returns tan _1 (6/a) e [— ir,ir], and, fur- 
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FIG. 3: Same as in Fig. Q but for the real damping. 



FIG. 4: Transition frequencies uji (rad/s) and u>2 (rad/s) for 
Rio = 1 lira, R20 = 4 /j,m, and the real damping, normalized 
by W10 (rad/s). 



thcrmore, adopted the operation 

j _ / atan2(A,-,B.,) + 2tt if atan2(A,-, Bj) < 0, 
i ) a,t&n2(Aj, Bj) otherwise 

in order to obtain results for Z ^ 1 that are consistent 
with the established knowledge of single-bubble dynam- 
ics, e.g., 4>j w 7T when u> > ujjo and Sj ss (j = 1 or 
2). 

The dashed curves displayed in Fig. Gfb) show results 
for w = 1.03^20 (= 0.413o;io), i.e., for ujiq > u> > u>2o- 
In this case, we can observe only one sharp shift of (j>\ 
at I w 12, causing the sign reversal. This shift of 4>\ is 
due to the second-highest transition frequency of bubble 
1 (this frequency also not corresponding to the natural 
frequency!), because the lowest transition frequencies of 
both the bubbles decrease as I decreases. 

These results reveal that in the above cases, the transi- 
tion frequencies other than the natural frequencies cause 
the sign reversal of the secondary Bjerknes force. This 
conclusion is obviously different from the previous inter- 
pretations described by means of the natural frequencies 

fiiiam. 

It is interesting to point out that, in the case where 
uj > cjio > uj20 and uj « luiq, the phase delay of the larger 
bubble was sometimes greater than it (see Fig. Ela))- 
Such a result cannot be given by a single-bubble model 
which predicts a phase delay of up to 7r. This may be 
explained as follows: When uj > uj\q > LO20 is true and 
I is sufficiently large, both bubbles pulsate out-of-phase 
with p cx , emitting sound waves whose phases are also 
out-of-phase with p CK . As / decreases, if uj w ujiq, the 
amplitude of the sound wave emitted by bubble 1 mea- 
sured at r 2 can be greater than the amplitude of p ox . In 
this situation, bubble 2 is driven by a sound wave whose 
oscillation phase is delayed by almost n from that of p cx . 
This results in 02 > 7r, because the pulsation phase of 
bubble 2 delays further from that of the sound wave. 

We show the here results given by using Eq. (|l(Jfl in 
order to examine the influences of the damping effects on 
the sign reversal and phase shifts. Figure shows the 
recalculated transition frequencies. As already discussed 
when the damping effects are not negligible, the 
bubbles have only one transition frequency in the large- 
l region. The solid curves displayed in Fig. [3 show <pi , 



(a) to = 1 .01 e> 10 (b) m = 1 .03 <a 
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/ ; 

FIG. 5: 4>i/tt, 02/t, and cos(<^i — (f>2) for (a) uj = l.Olwio 
(rad/s) and (b) uj = 1.03oj20 (rad/s). The dashed curves 
denote the results for the reduced damping. 



4>2, and cos(</)i — fo) for uj — l.Olwio and uj — 1.03^20- 
Their tendencies are similar to those given with the re- 
duced damping, although their profiles are smoothed sig- 
nificantly (Such a smoothing of the phase change by the 
damping effects is well known for a single-bubble case) 
and the points at which the sign reversal takes place are 
shifted slightly; the positions of these points are, in the 
case of uj = l.Ol^io, I ~ 13.59 for Sj and I w 14.57 for 
Sj/100, and, in the case of uj — 1.03^20, I ~ 12.66 for Sj 
and I w 12.48 for 5j/100. Moreover, 4> 2 for uj = LOl^io 
does not exceed 37r/2 (the minimum value of UJ2 larger 
than uj\o is 1.027wio.); even so, the sign reversal occurs at 
almost the same point as that given with <5j/100, away 
from the point where <p\ — tt/2. This result may be 
interpreted as the "vestige" of the highest transition fre- 
quency of the larger bubble having given rise to this sign 
reversal. Detailed theoretical discussions for the slight 
shift in / for cos(0i — <\>2) = due to the damping effects 
will be provided in a future paper. 
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(a) Natural frequencies [D-Z's assumption] 
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FIG. 6: Comparison between the theoretical and the numer- 
ical results. The lines and the circles denote the theoreti- 
cal and the numerical results, respectively, of cos(0i — (j>2) 
for uj — l.Olwio (rad/s) and 1.03W20 (rad/s). Additionally, 
sgn^ViVb)) for P a /Po = 0.2 (the dots) and 0.5 (the dash- 
dotted curves) are plotted for a brief investigation of nonlinear 
effects. 



Next, we show results for smaller bubbles (Rio — 1 /mi 
and i?2o = 4 ^m). The value for viscous loss is used for 
the damping coefficients, i.e., 



4p 



(ii) 



where the viscosity of water p, = 1.002 x 10~ 3 kg/(m s). 
Because the thermal effect is neglected, n = 7 = 1.4. 
Figure 0] shows the transition frequencies, and Figure 
shows </>i, 02, and cos(0i — §2) for to = I.OIojio and 
uj = 1.03w 20 (= 0.201wio) with <5,/100 (the dashed 
curves) and Sj (the solid curves). The qualitative natures 
of these results are quite similar with the previous ones; 
thus, additional discussion may not be necessary. Using 
this example, we perform here a comparative study of the 
theoretical results with the numerical results in order to 
confirm the former's correctness. In the numerical experi- 
ment, we employ the coupled RPNNP (Rayleigh, Plesset, 
Noltingk, Neppiras, and Poritsky) equations (see, e.g., 
Ref. E3); 



" 3 - 2 1 

ivL-Rl + -itj p w .i — 

2 p 

R2R2 + TT-Rf Pw,2 = 

2 p 



1 


Pox 4 


p d 


p 


' Ddt 


1 


PCX ' 


p d 


p 


^ Ddi 



where 



P 



w,0 



Pn 



2a 
Rjo 



Rjo_ 
R, 



3 k 



2(7 

^7 



This system of nonlinear differential equations are solved 
numerically through the use of the fourth-order Runge- 
Kutta method in which i?2, Ri, and R2 are used 

as dependent variables, and (r\R\R\R2^ [oc (ViV^j 

in Eq. 10] is then calculated. The time average is 
performed during a sufficiently large period after the 

transients have decayed. Normalizing (RIR1R2R2 



(b) Natural frequencies [Theory] 
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(c) Transition frequencies [Theory] 
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FIG. 7: Characteristic frequencies of two coupled bubbles and 
different interpretations of the sign reversal. The dashed lines 
show a typical driving frequency lying between wio and u>20, 
where uiio > u>20 is assumed. Doinikov and Zavtrak assumed 
that the natural frequencies of both the bubbles increase as 
the bubbles approach each other (a). Assuming this, the sign 
reversal for uiio > u > CJ20 seems to be explained. This 
assumption is, however, inconsistent with the classical theory 
for natural frequencies (b). The present theory can explain 
this reversal without such an inconsistency (c). 



yields the numerical approximation of cos(0i — (^2), 
where max(\Rj(t) — Rjo\) indicates the pulsation ampli- 
tude of bubble j given numerically. The amplitude of the 
external sound is set to P a — O.OlPo- In Figure|51 the nu- 
merical and the theoretical results are displayed in piles. 
These results are in excellent agreement, confirming the 
correctness of the theoretical results given above. In the 

same figure, we have shown additionally sgn^ViV2^ 

for P a = 0.2P (the dots) and 0.5P (the dash-dotted 
curves) in order to briefly investigate nonlinear effects 
on the sign reversal, where sgn(A") = 1 for X > and 
sgn(AT) = —1 otherwise. In plotting these results, we 
omitted the data in the case where i?i(i) + R2(t) > D 
was observed during the computation. As is clearly 
shown, increasing the driving pressure reduces the 
distance for which the sign reversal takes place. This 
result appears to be consistent with the well-known 
nonlinear phenomenon that a strong driving pressure 
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decreases a bubble's (effective) resonance frequency, see, 
e.g., Refs. @,0| . (Imagine that the transition frequencies 
shown, e.g., in Fig. decrease but the driving frequency 
holds, which might shorten the distance for F = 0.) 
More detailed and concrete discussions on the nonlinear 
effects will be provided in a future paper. 

To summarize our discussion, we compare the present 
interpretation with the previous ones. FigureEfa) shows 
the dependency of natural frequencies on I that Doinikov 
and Zavtrak assumed 0, Il0l | . Their assumption explains 
the sign reversal occurring when u)±o > to > lu 2 q, f° r 
example, as taking place around I at which lu 2 = lu is 
true. Yet, as mentioned, their assumption is inconsistent 
with the theoretical results regarding natural frequencies 
given previously 0,0 (see Fig.0[b)). On the other hand, 
it is difficult to determine by only observing the natural 
frequencies that the sign reversal can take place for luiq > 
lu > W20, because the classical theory does not show that 
a kind of characteristic frequency exists in the frequency 
region between luiq and lu 2 q ■ The present theory explains 
the sign reversal in this case as taking place around I at 
which lui — lu is true (see Fig. 0c), where we assume for 
simplicity that the damping effect is negligible), and is 
consistent with the theory for natural frequencies because 
the transition frequencies include the natural frequencies. 

IV. CONCLUSION 

We have investigated the influences of change in the 
transition frequencies of gas bubbles, resulting from their 
radiative interaction, on the sign of the secondary Bjerk- 
nes force. The most important point suggested in this 
paper is that the transition frequencies that cannot be 
derived by the natural-frequency analysis cause the sign 



reversal in the cases of both ui > W10 > lu 2 o and 
lu w > lu > lu 2 q- This interpretation has not been pro- 
posed previously. The present results also show that the 
theory given in Ref. |l6| for evaluating the transition fre- 
quencies of interacting bubbles is a reasonable tool for 
accurately understanding the mechanism of this reversal. 
In a paper currently in preparation p3. w e will use the 
direct numerical simulation technique |2J, |25| to verify 
the present theoretical results. 

Lastly, we make further remarks regarding the results 
described in Ref. ^(|- I n that paper, the frequency of the 
external sound (/ = lu/2-k) was assumed to be / = 63 
kHz, which is 60 times higher than the partial resonance 
frequency of a bubble of Ro = 3 mm (1.094 kHz); never- 
theless, the reversal was observed at a very small I. (In 
Ref. |9j , the driving frequency is assumed to be compara- 
ble to the partial natural frequencies of bubbles, and the 
bubble radii are several tens of micrometers.) This result 
reveals implicitly that the mathematical model proposed 
in Ref. [3, which takes into account the shape deviation 
of the bubbles, predicts such a strong increase of the tran- 
sition frequencies of closely coupled large bubbles that 
this increase cannot be explained by the classical model 
for coupled oscillators used here. Derivation of the tran- 
sition frequencies of Doinikov and Zavtrak's model would 
be an interesting subject for future study. 
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